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Abstract 

A Dirac picture perturbation theory is developed for the time evolution 
operator in classical dynamics in the spirit of the Schwinger-Feynman-Dyson 
perturbation expansion and detailed rules are derived for computations. 
Complexification formalisms are given for the time evolution operator suitable 
for phase space analyses, and then extended to a two-dimensional setting 
for a study of the geometrical Berry phase as an example. Finally a direct 
integration of Hamilton's equations is shown to lead naturally to a path 
integral expression, as a resolution of the identity, as applied to arbitrary 
functions of generalized coordinates and momenta. 
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1 Introduction 



E. B. Manoukian and N. Yongram 



With the rapid progress of field theory in describing the basic interactions occurring 
in nature, it is becoming more and more important to extend the powerful techniques 
of field theory in describing physics, in general, in a more unified language. Such a 
unification has always been the goal of physics research over the years. The present 
work is aimed in describing classical dynamics in the spirit of field theory methods. 
Some earlier papers on this subject are given in Abrikosov (1993) pQ, Gozzi et al. 
(1989) 0, Schwartz (1976) @], and Wetterich (1997) [6 , which are, however, only 
tangentially related to our present investigations. We develop a Dirac (interac- 
tion) picture perturbation theory in the spirit of the Schwinger-Feynman-Dyson 
perturbation theory to all orders in a coupling parameter and detailed rules are 
derived for computations starting from the time evolution of any function of the 
generalized coordinates and momenta. Complexifications of the time evolution are 
developed suitable for phase space analyses and then extended to a two-dimensional 
setting to describe the so-called geometrical Berry phase (Berry, 1984 2\; Shapere 
and Wilczek, 1989 5 ) as an example. Finally, since the classical limit of the path 
integral, starting from the quantum regime, reduces to just a phase factor involving 
the classical action, such a limit (Abrikosov, 1993 [1 ; Gozzi et al, 1989 |3|) is 
not of very practical value, as it stands, for actual computations. Instead, we 
develop a path integral expression by direct integration of Hamilton's equations, as 
a resolution of the identity, that may be applied to any function of the generalized 
coordinates and momenta, in the same spirit of developing the resolution of the 
identity of a self-adjoint operator, in quantum physics, that may be applied to 
any vector in the underlying Hilbert space. The perturbation expansion, with 
the derived rules, is given in Section El while the complexification formalisms are 
developed in Section El Section |U deals with the path integral expression and its 
consistency with the Poisson-bracket solution. 
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2 Rules for Computations and the Dirac Picture 
Perturbation Theory 

The time derivative of an arbitrary function f[q(t),p(t)] of generalized coordinates 
and generalized momenta, via the Poisson-bracket formalism, reads 



-f[q(t),p(t)} = 



dH(t) d dH(t) d 



dp(t) dq(t) dq(t) dp(t) 
fMt),p(t)] 



f[q(t),p(t)} 



(2.1) 



where H[q(t), p(t)] = H(t) is the Hamiltonian, constructed out of the variables q(t) 
and pit), assumed with no explicit time dependence. Similarly, 



f[q(t),p(t)]=fn[q(t),p(t)] 



dH(t) d dH(t) d 



_ dp(t) dq(t) dq(t) dp(t) 
leading to the familiar time evolution 

f[q(t),p(t)}=e td f[q,p} 

where 



O 



dH d dH d 



dp dq dq dp 



f[q(t),p(t)} (2.2) 

(2.3) 
(2.4) 



and q = q(0), p = p(0), H = H(0). 
Upon using the integral identity 

e t[A + XB] e -tA = 1 + X [\ Te T[UxB]z e -TA 

Jo 

for two operators A and B, and setting 

H = H 1 + \H 2 

'dHx d dH 1 d 



A 
B 



dp dq dq dp 
dH 2 d _ dH 2 d 
dp dq dq dp 



(2.5) 

(2.6) 
(2.7) 

(2.8) 
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one obtains, upon iteration of (j2.5j) . the expression 

ft ft\ ftn-l 



f[q(t),p(t)} = J2 Xn dtl di 2---/ <tt n B{t n )---B{t x )e tA f[q,p] (2.9) 

n=0 - 70 70 J ° 



where 

5(t) = e tA Be~ tA 
with -B(i) independent of the parameter A. 

In particular, for 



fl(f) = F 



q + —p 

m 



d t d 
dp m dq 



F[q] = -V'(q) 

ffap] = q 



then 



and we obtain 



e tA f[q,p] =q + —p 
m 



q(t) = V \ n / dtj / dt 2 ■ ■ • / dt n B(t n ) • ■ • Bfa) 

n=0 J ° J ° J ° 



q H p 

m 



(2.10) 



(2.11) 
(2.12) 

(2.13) 
(2.14) 

(2.15) 
(2.16) 



In detail 

B(* 2 )fl(fi) 
and so on. Upon letting 





t 




9 + 


— p 


= ( 




m 






* 




9 + 


— p 


= ( 




m 





i - 1 



m 

h-t 2 



<7 H p 

m 



/;/ 



t-t 



, *2 

g H p 

m 



F' 



9 H P 

m 



q H P = Ui 

m 



(2.17) 
(2.18) 

(2.19) 
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we obtain, after a systematic analysis, the following general explicit rule : 



B{t n )---B{t x ) 
't-ti 



t 

q H p 

m 



m 



x 



E 



F[u r 



m 



t 2 -h 



m 



<5(fci,2) 



<5(fc 2 ,3) 



t\ — t Tl 



m 



<5(fci ,n) 



^2 — t n 



m 



x 



^n— 1 tr 



III 



6(k n -l,ri) 



where 



and the sum in ()2.2())1 is over all fc's and 8's such that 

h + h + h k n -i —n — 1 

k\ = 1, . . . , n — 1 
fc 2 = 0,l,...,n-2 ' 



fc n _i = 0, 1 



and 



5(fci,j)=0, if k i = 
6(ki,j)=0, if Ui<i 

and for j > i, the 8(ki,j) are zero or one such that 



E j) = 

j=i+i 



i — 1, . . . , n — 1 



and 



71-1 



5]%J) = 1, (5(^,J)=0, j = 1,2,...,*) 

8=1 



(2.20) 



(2.21) 



(2.22) 



(2.23) 
(2.24) 



(2.25) 



(2.26) 
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that is, for a fixed j , tj appears only once in the product : 



3=2 3=3 



(2.27) 



For example, for n — 4, 

ki = 1, 2, 3; fc 2 = 0,l,2; fc 3 = 0, 1; fci + k 2 + k 3 = 3 

and for k\ = 2, k 2 = 1, A; 3 = 0, 

5(2, 2) + 5(2, 3) + 5(2, 4) = 2 
5(1, 3) + 5(1, 4) = 1 

and 

5(2,j)+5(1,j) = 1- 
For the harmonic oscillator V(q) = q 2 /2, A = muj 2 , and 



~ ^ d ^ v d 



m dq 



The only solution being 
with 



k 1 = ... = k n _ x = 1 
F^[ Ul ] = -l, i = l,...,n-l; %] = - 



g H p 

m 



we obtain 



?(*) 



00 rt rt\ rtn-i 

n=Q Jo Jo Jo 



X (ti — t 2 ) ■ ■ ■ (t n -i — t n 

which integrates out to 

oo 

q(t)=qJ2(-" 2 ) 



q H p 

m 



j.2n 00 
t M p 



n=0 



(2n)\ m 



n=0 



(2n+ 1)! 



or 



q(t) = q cos out + 



p sin cut 



(2.28) 

(2.29) 
(2.30) 

(2.31) 

(2.32) 

(2.33) 
(2.34) 



m uj 



(2.35) 

(2.36) 
(2.37) 
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3 Complexification of the Time Evolution 

We consider the complex dynamical variable 

Z{t) = aq{t) + iftp(t) 



(3.1] 



where a and b are arbitrary real constants. Equation ()2.9j) then immediately leads 
to 



00 r-t rti 

Z(t) =^A n / dt x / dt 2 



n=0 



JO 



dt n B{t n )B{t n ^) 



x B(t! 



t 



a q-\ p + ibp 



m 



(3.2) 



in which B{t n ) ■ ■ ■ B{t\) [a(q + tp/m) + \bp\ coincides with the expression in (|2.2 
except that the overall factor (t — t\)/m in the latter is replaced by 



a 1 ) + 10 

m 



(3.3) 



A more explicit expression may be also obtained directly from ()2.3)1 . 
end, consider the Hamiltonian 



2 2 

H = ^- + mu% + \V{q) 
2m 2 



and upon setting 
one readily obtains 



Z(t) = q(t) + i 



:P(t) 



muo 



Z{t) = exp 



2 



Z(0) 



where 



d d 
C=(Z-Z*)[^ + 



dZ dZ* 

2A -UZ + Z* 



To this 
(3.4) 

(3.5) 
(3.6) 



d d 
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where F[. . .] is defined in fl2.13|) . For V(q) = —q, 
C ■ 



muj z I oZ \ muj z I oZ* 



leading from 



Z(t) = exp 
and upon using the identity 

exp 



-itu) I Z 



7 9 

aZ— — 

dZ 



for any constant a, to the expression 

Z{t) x 



A \ d 



mui 2 ) dZ 



Z = Z exp a 



e~ ltuJ Z 



muo A 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



A more interesting application is to the geometrical phase associated with the 
famous Foucault pendulum with Hamiltonian 



(3.12) 



where \fgjL = lu ^> lu, u z = cjsinA, with A denoting the latitude and L the length 
of the pendulum. Then 



d =i^ + y^)^-( m ^ x -py^ dPx 

fPv \ 9 , 2 \ 9 

+ Vm ~ XUJz ) dy ~ \ muJ ° y + PxUJz > 







d P y 



(3.13) 



Upon defining 



U 
V 



x 



X + 



rn 

iPx 
rn 



> I y + 
1 1 y + 



iPy 

m 
m 



(3.14) 
(3.15) 
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then 



and 



Z = x + iy 



= 



1 



U+ V* 



U + UJ Z 



(3.16) 



XT 



d 



dU* 



i dV 



dV* 



which lead to 
Z(t) = exp 
With initial conditions 



it (u + U z )U-^j 



exp 



it (cg> - oo z ) V* 



d 



dV* 



U + V* 



x(0) = e^x , 2/(0) = 0, x(0) = 0, 2/(0) = 



(3.17) 



(3.18) 



(3.19) 



where e l2n denotes the initial phase of the plane of oscillations of the pendulum 
moving clockwise, we have 



1/(0) = e l27T x = V*(0) 



giving 



Z(t) = e- itu)z e i2n x cos to t. 
For t = 2ir/u, corresponding to a full rotation of the earth, 



Z\ — ] = exp 



i27r (l — sin A) 



Xq cos Ztt — 



(3.20) 
(3.21) 

(3.22) 



giving the familiar geometrical Berry phase exp [i27r(l — sin A)], which is indepen- 
dent of g and u, with 27r(l — sin A) denoting the solid angle subtended at the center 
of the earth. 



4 A Path Integral Expression 

Upon integrating Hamilton's equations dH/dp = q, dH/dq = —p we have 



f[q(t),P(t)] = f 



,(0)+fd,fW 
Jo dp{V 



p(0) 



dq{V 



(4.1) 
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We devide the time interval from to t into n subintervals each of length t/n. 
Evaluating each of the summands at the left-end points of these subintervals, we 
may rewrite the right-hand side of ()4.1j) as 



lim / 



t dH k t dH k 

Q0 n{^dp k ' P0 n{^ Q dq k 



(4.2) 



where now we set q(0) = qo,p{Q) = Po, and & is a short-hand for tk/n. The expression 
(14. 2J) to which the limit n —>■ oo is to be taken may be equivalently rewritten as 

/ dq n dp n 5 ( g + - ^ W~ Hk ~ qn ) S [ Po ~ ~ 7f~ Hk ~ Pn I ^ gn ' Pn ^ 
J \ n k =o ° Pk / \ n k=o ° qk J 



\\ dq k dp k 5 I q k 



k=l 



k=0 

c-l i — — — qk 



ndp k - 



i 



si 1 9 TJ 

X(5 p k -i iifc_i - p k 



n dq k -i 



j J dq k dp k 5 



k=l 



t 



d 



n { dp 



■Hk-i 



f[ln,P, 
(q k - q k -i) 



k-i 



t/n 



x,)[ - 



n [ dq 



d tt , (Pfc ~ Pfc-1; 



Jfe-1 



t/n 



f[Qn,Pn}- 



(4.3) 



Upon taking the limit n — > oo of (|4.3)1 . Eq. (|4.1jl becomes 

/fe(f),p(f)] = J 9[q) 9\p] 5 {^-i) 5 {^ +p ) ( 4 - 4 ) 

as a resolution of the identity as applied to arbitrary functions f[q(t),p(t)] of g(t) 
and p(t). The integrals are over all paths in phase space, from to t, starting 
from (g(0),p(0)), and the delta functionals as obtained directly from Hamilton's 
equations restrict these paths to the classic one obeying Hamilton's equations at 
each instant of time. Finally, we have also used the fact that q n — > q(t), p n — > p(t) 
from the very definitions in (|4.2|) . ()4.3|) for n — ► oo. Equation (|4.4jl is in the spirit 
of the resolution of the identity of a self-adjoint operator in quantum physics as 
applied to any arbitrary vector in the underlying Hilbert space. 
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For applications, (J4.4j) may be rewritten in the more manageable form 



9[q\ 9\p\ ®[\) @[rj[ exp 



+t)(t 



dp(r) 
dH(r 



T 



+ p(r] 



f[q(t),p(t)) (4.5) 



in terms an uncountable infinite number of Lagrange multipliers A(-), //(•). 

For example, for a charged particle of charge e in a uniform magnetic field B, 
say, along the z-axis [B = (0,0, B) ], we may write for the vector potential A = 
( — <?2, qi) B/2, with motion in a plane, and for the Hamiltonian H 



TT 1 / mui\ 2 1 / mui\ 7 



(4.6) 



with uj = eB/mc, A = (Ai, A2), r] = (771, 772)- The time-integrand in the exponential 
in 1)4.51) . without the i factor, may be rewritten as 



Ai 

+ A 2 



1 / muj 
— Pi + 12^ 



m 



Qi 



+ Vi 



lu ( into 



1 / muj 

— \P2-q1 — 
m V 2 



A x 
2mu 

m 



Q2 



+ V2 



uj ( muj 



+ Pi 

+ P2 



+ 



4i 



dr 
d 



dr 



+ iu)U 



dr 



iu)U* + 4-(V - V 
dr 



A? 



2muj\ 



where 



02 

4 



U 
V 



dr 



dr 



iu U 



dr 



h)\U* - —{V-V* 
dr 



dr 



-\uo\U* - — (V + V 
dr 



mw 



- P2 + 1 [Pi + <?2 



2 r V V " 2 
ma; \ / rmo 
Qi-pr +P2 1 + 1 [Pi ~ Q2-^r- 



(4.7) 

(4.8) 
(4.9) 
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The coefficients of Ai, 771, A2, ^2 on the right-hand sides of (|4.7|) are all reals. Upon 
integration over Ai, rji, A 2 , t] 2 , we learn that the real and imaginary parts of 



d ■ \ tt , d Tr 
— + iu )U± —V 
dt J dt 



must vanish. That is, V(t) = V(0), and (d/dt + iw) U = or U(t) = U(0)e~ 
Upon equating the real and imaginary parts of each of the latter two equations we 
obtain the solution 



92 (*) 



gi(Q) , P2(0) ) 
2 mu 

'g 2 (0) + pi(0) 



9i (0) pa(0) 



2 mcj 
, sin tut 

2 ma; 
92(0) pi(0)\ , /g 2 (0) , Pl (0) 



cos cut 



+ 



2 mcu / \ 2 

qi (o) , P2(o) 



(4.10) 



cos cut 



+ 



sin ujt. 



(4.11) 



To check the consistency of (|4.4|) . we note that in reference to the first equality 
on the right-hand side of (|4.3|) . that 



cxp 



t_fdH k ^\ 

n 



dpk-i J dq, 



k-i 



is not quite a translation operator for a function of qk-i, since dHk-i/dpk-i may, in 
general, depend on qk-\ as well. However, in view of the fact that the limit n — > 00 
is to be taken, this operator may be indeed taken to have such a property for the 
accuracy needed. A similar comment applies to the 



exp 



t 

n 



ggfc-i 



5 



fc-i 
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operator. To the accuracy needed (|4.1jl - (|4.Hjl lead to 

t fdH k -x\ d 



f[q(t),p(t)] = Hm Jy Y[dq k dp k jexp 



x exp 



n V dpk-i J dq k -i 
t (dRk-i 



d 



xS{q k -i -qk)S(Pk-i -Pk)\ I f[q n ,Pn\- (4.12) 



Upon integration over (q k ,p k ) with the aid of the delta functions which eventually 
pick up the (qo,Po) value for the former, we obtain the Poisson-bracket solution 



f[q(t) ,p(t)] = lim exp 





~tdH(0) d ' 




exp 


n dp(0) dqo 


exp 



exp 



f dH{0) d _ dH{0) d 
\ dp(0) dq dq(0) dp 



tdH(0) d 
n dq(0) dpo 

f[qo,Po]- 



f[qo,Po] 



(4-13) 
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